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1. Introduction 

The invention of a transfer matrix technique for dynamical triangulations [II1J2] has 
made it possible to introduce a concept of distance in two-dimensional quantum gravity. 
This again has made it possible, in the case of pure gravity, to calculate analytically a loop- 
loop correlation function as a function of geodesic distance and to show that the scaling 
relations from flat space statistical mechanics have analogues in the quantum theory || . An 
intriguing result of these investigations is that the fractal dimension of the quantum space 
times turns out to be four and not two. Unfortunately, when it comes to the inclusion of 
matter fields the transfer matrix technique for dynamical triangulations has not yet proven 
as efficient as its regular lattice counterpart. Correlators as a function of geodesic distance 
have been studied numerically 0-0] but analytical calculations have not been possible.^ 




> 
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Fig. 1: A typical discretized universe of (l + l)-dimensional Lorentzian quan- 
tum gravity. Each constant-time strip is an arbitrary succession of triangles 
pointing up or down. Boundary conditions can be chosen to be periodic, free, 
fixed according to the model at hand. 

Recently, a different discrete approach to two-dimensional quantum gravity has been 
proposed |1^,|13]]. In this approach which is known as Lorentzian quantum gravity only 
triangulations which admit a causal structure are allowed in the state sum. An example of 
such a triangulation is shown in Fig. 1. The elementary building blocks are triangles with 
one space-like and two time-like edges which are glued together at random so as to form 



4 There have, however, been such calculations using an alternative, fully valid, but not equiv- 
alent definition of distance |8-1C], see also [11]. 
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a piling of constant-time slices. In a sense such triangulations represent only the space 
fluctuations of the universe while respecting chronology. In Lorentzian quantum gravity, 
baby universes are not allowed. In other words, a given time slice of the triangulation 
has only one connected component. This, in particular, means that a universe always 
has spherical topology. The model of pure Lorentzian quantum gravity was solved in ||12|| 
exploiting the fact that the generating functional for its transfer matrix obeyed a certain 
iterative equation. The continuum limit of the model proved to be different from that of 
usual Euclidean quantum gravity or Liouville theory. For instance, the fractal dimension 
of the quantum universes turned out to be equal to two. However, the continuum limit 
of Lorentzian quantum gravity is in accordance with continuum results obtained from 
calculations carried out in the proper time gauge |14[ .1 A Lorentzian lattice by construction 
is something between a regular lattice and a truly random one. It is therefore interesting to 
know if the introduction of matter fields on such lattices leads to any interaction between 



matter and geometry. In |15[ the effect of coupling an Ising model to Lorentzian quantum 
gravity was studied. The investigations were based partly on a high temperature expansion, 
partly on Monte Carlo simulations. No sign of interaction between matter and geometry 
was seen. Numerical simulations currently being carried out indicate, however, that for 
matter fields with c > 1 a non-trivial interaction between matter and geometry takes place 

In a first part of the present article we formulate and solve exactly various other 
extensions of two-dimensional Lorentzian quantum gravity. The simplest extension, which 
we develop in details throughout the paper, is that of Lorentzian gravity with a higher 
curvature term which is equivalent to Lorentzian quantum gravity coupled to a simple 
dimer model or to Lorentzian quantum gravity with universes built from triangles and 
squares. We show that this model is integrable in the sense that its transfer matrix can be 
explicitly diagonalized and that there exists a set of commuting transfer matrices which 
allow us to interpolate between regular and random Lorentzian lattices. Furthermore we 
describe how to construct models of Lorentzian quantum gravity allowing general p-gons 
as building blocks while preserving the integrability structure. As will become evident 
the transfer matrix formulation is by far more efficient for Lorentzian than for Euclidean 
triangulations. In fact the integrability structure revealed has many treats in common 



If baby universes are allowed in Lorentzian quantum gravity the usual Euclidean gravity 



continuum limit is obtained. We refer to [12] for a discussion of this point. 
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with the integrability structure found in regular lattice transfer matrix studies and it is 
our hope that our further investigation of this structure will enable us to solve exactly more 
realistic systems of matter fields coupled to Lorentzian quantum gravity such as Lorentzian 
triangulations equipped with Ising spins. For all the models that we have considered so 
far the universality class of the geometrical system is the same as that of pure Lorentzian 
quantum gravity 

In a second part of this article we study the relationship between Lorentzian quan- 
tum gravity and random walks. We prove that there exists a one-to-one correspondence 
between Lorentzian triangulations and directed random walks drawn on the regular tri- 
angular lattice. This allows us to set up a dictionary connecting concepts in Lorentzian 
quantum gravity to concepts in the theory of random walks. For instance, the expression 
for the loop-loop correlation function follows from a similar expression for large excursion 
probabilities for random walks. Furthermore, the integrability structure of the model of 
Lorentzian quantum gravity with a higher curvature term can be understood in terms of 
the possibility of a simple block decomposition of the corresponding random walk. Finally 
the random walk equivalence provides an explanation why Lorentzian triangulations have 
fractal dimension two and why we can not obtain from our model with a higher curvature 
term other critical behaviour than that of pure Lorentzian quantum gravityi 

The simple extension of Lorentzian quantum gravity involving a higher curvature term 
as well as its various alternative interpretations are described in Sect. 2. Sect. 3 treats 
the equivalence between random walks and Lorentzian triangulations and Sect. 4 contains 
a discussion of Lorentzian gravity involving general polygonal building blocks. Finally, in 
Sect. 5 we conclude and discuss the future prospects of transfer matrix techniques for 
Lorentzian triangulations. 



2. Discrete Lorentzian 2D Gravity via Triangulations 

2.1. Lorentzian Quantum Gravity with a Higher Curvature Term 

We consider quantum universes of the type depicted in Fig. 1. Generally we will be 
interested in a time evolution from, say, time t = 1 to t = T, i.e. a time strip of width T. 



6 Dimensional analysis indicates that adding a higher curvature term to the Einstein-Hilbert 
action should not modify continuum physics. However, this argument is only valid perturbatively 
and one could still hope for the existence of a (non-perturbative) ultra-violet fixed point, a scenario 



(in 4D) denoted by Weinberg as asymptotic safety [17|. 
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In the course of the article we will consider various possible boundary conditions in space. 
These will be explained at the relevant points. 

As usual in two-dimensional gravity, we will attach a weight g per triangular face, 
resulting in an overall factor of g A for each triangulation, where A is the total number of 
triangles, measuring the area if we decide that triangles have unit area. 

We now generalize the model by also attaching to each triangulation an intrinsic cur- 
vature weight defined as follows. For a two-dimensional triangulated manifold, curvature 
resides on vertices, and the curvature at a vertex is proportional to (v — 6) where v is the 
valence of the vertex. Here we do not wish to add an ordinary curvature term (which is 
trivial in two dimensions) but a term which suppresses (or enforces) local curvature. Due 
to the time slice structure of the triangulation we take as a measure of the curvature at a 
vertex the quantity \v 1 — 3| + |t>2 — 3| where v± and i>2 are the number of triangles adjacent 
to the vertex in the upper and the lower time-slice respectively, and we attach to each 
vertex of the triangulation the weight a^ 1-3 '^ 2-3 '**/ 2 . Introducing this weight factor can 
be viewed as adding a higher curvature term to the Einstein-Hilbert action (note, however, 
that one can have v — 6 = and \v\ — 3| + \v2 — 3| ^ 0). Our model is similar in spirit to 



the models considered in |18|], where the effect of adding a higher curvature term to usual 
dynamically triangulated gravity was considered. Since any vertex is linked to exactly 
two up (resp. down) triangles is its upper (resp. lower) slice, with an arbitrary number of 
down (resp. up) triangles in-between, a simpler, completely equivalent way to introduce 
the curvature weight is to attach a weight a per pair of adjacent triangles pointing both up 
or both down, within the same constant-time slice. The graphical representation of Fig.l 
is more transparent in the dual picture, where a constant-time strip becomes a succession 
of half-edges attached to the dual constant-time line, where the half-edges lying above 
(resp. below) the line correspond to triangles pointing down (resp. up). Moreover, as each 
triangle has exactly one space-like edge, the triangles pointing down in the strip of time t 
all share an edge with the triangles pointing up in the strip of time t + 1, henceforth the 
dual half-edges are connected to form edges between the lines of time t and t + 1. 

The Boltzmann weights defined above are depicted in Fig. 2, where we also give an 
example of a world-sheet configuration in the dual picture. 

The Lorentzian gravity problem above, and its generalization with intrinsic curvature 
admit a transfer matrix formulation. Indeed, if we compute the Boltzmann weight of a 
single time-line with i lower half-edges and j upper ones, and sum over all possible relative 
positions of these half-edges, we find a number Tij(g, a) that depends only on the numbers 
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Fig. 2: The Boltzmann weights of the Lorentzian gravity model with a 
higher curvature term. Each vertex (of the dual lattice) receives a factor g. 
In addition each sequence of neighbouring up-up or down-down half-edges 
receives a factor a. We have represented a typical world-sheet configuration. 

of half-edges, but not on their specific relative positions. This can in turn be taken as the 
transfer matrix element that transfers from a row of i half-edges to one of j. For 

instance, the partition function of a strip of time width T, with l\ lower triangles pointing 
up and £2 upper ones pointing down is simply 

Zz lji2 (T\g,a) = (T(g,a) T ) . (2.1) 

V / e 1} £ 2 

For periodic boundary conditions, this is the so-called two-loop correlator of Lorentzian 
gravity, which describes world sheets of cylindric form bounded by two loops of respective 
lengths £\ and £2- 

Let us adopt the following choice of fixed boundary conditions, that differ from the 
periodic ones of |L3|]. We assume that each time strip always has one lower half-edge 
(i.e. triangle pointing up) on its leftmost end and one upper one (triangle pointing down) 
on its rightmost end. This implies in particular that i, j > 1 in the transfer matrix, i.e. 
that the time slices never degenerate into the vacuum. These boundary conditions can be 
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described as "staircase" conditions since they imply that the left and right hand side space- 
like boundaries of the world-sheet have the shape of a staircase with the stairs extending 
to the right at each successive time line (see the example of Fig. 2). For reasons that will 
become clear later, we attach a weight ^fg only to each of the two boundary half-edges in 
each slice. 

We may now compute the transfer matrix element Tij(g,a) by summing over all 
possible configurations of upper and lower half-edges. This is done by first summing over 
the number k > 1 of blocks made, say, of n r > 1 consecutive lower half-edges, followed by 
m r > 1 upper ones, r = 1, 2, /c, and then summing over the corresponding partitions of 
i = J2 r n r and j = J2 r m r-> with positive integers n r , m r . Note that the conditions ni > 1 
and rrik > 1 are a direct consequence of the boundary conditions we have imposed. We get 



mi 



m 2 



fc^l n r ,m r >l , r=l ,2 , . . . ,k 



n-2 



= E E 



^i+j-l a S(n r -l)+E(m 7 .-l) 



/c^l n r ■, 1Tl r ^ 1 ■, T = l,2,...,fc 

= n t+ i- 1 a i+J V — 
k>l v 



3-1 



(2.2) 

This is in general a polynomial of a, as the range of the summation is finite (k < min(i, j)). 
The formula (|2.2| ) is best encapsulated into the generating function 

0{x,y\g,a) = ^ x'y J Tij(g : a) 

_ 1 ^ 1 

(gax) k (gay) k 



fc>l i>fc V 7 j>/c V 



-Y-- 

g ^ a 2k (1 



(2.3) 



(1 - gax) k (l - gay) 1 



gxy 



1 — ga(x + y) — g 2 {l — a 2 )xy 
In the particular case a = 1, this reduces to the generating function for pure Lorentzian 



gravity discretized via triangulations to be compared with that of [[121 , but with different 
boundary conditions. When a = 0, ( |2~3| ) reduces to the generating function for flat world- 
sheets. Indeed, 9(x,y\g, 0) = gxy /{I — g 2 xy) simply means that Tij(g, 0) = g' l+ i~ 1 5ij. 
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In other words, the pattern of triangles becomes regular (successions of triangles pointing 
up,down,up,...etc...), and the discrete world-sheet becomes a rhombus-shaped portion of 
the triangular lattice, with area 2£ x T, where £ denotes the number of lower triangles 
pointing up in the first time slice. Hence our parameter a allows us to interpolate between 
flat and curved spaces. 



2.2. Diagonalization of the transfer matrix 

To solve our model, and be able to compute quantities like (|2.1|) , the simplest strategy 



is to diagonalize the transfer matrix T(g, a). But as the number of triangles in each time 
slice is not bounded, the matrix has infinite size. Nevertheless, at the expense of extra 
precautions, we have been able to construct a complete set of orthonormal eigenvectors for 
T(g, a), as a linear map acting on the infinite-dimensional vector space E = JR ® 1R ® .... 
They read as follows. 

Unless otherwise noted, we will always assume that the parameters g and a are real 
and satisfy the inequalities 

,l-^ 2 (l-a 2 ), , , , . 

^ q>2, ga<l. 2.4 

ga 

Let q = q(g,a) be the real solution to the quadratic equation 

ga(q+-) = l-g 2 (l-a 2 ), (2.5) 
Q 

such that |g| < 1. Then the functions 

F m (x\q) = J2^ m) (Q) 



i>i 

\m—l 



,2 



x(q — x) 



(2.6) 



'1 — qx) 1 



for m = 1, 2, 3, ... are the generating functions for the components of the m-th eigen- 
vector = v( m '(g, a) for T(g, a). Moreover the corresponding vectors form an orthonor- 
mal basis of E. 



The second statement is readily proved by considering the following contour integral 
over the unit circle 

v (m+p) . v (m) _ V^(m+p)^(m) 



£ 

i>l 

f J^L Fm + ( x \q)F m (-\q) 
J ztirx x 



= -(l-q 



dx (q — x) p 1 ^ 



2z7r (1 - 

if p > 1 

1 if p = 

where, using the Cauchy residue formula, we have noted that there was no pole of the 
integrand inside the unit disc when p > 1 (recall that \q\ < 1), and the result for p = is 
simply given by the residue at the pole x = q. This proves the orthonormality of the set 

(m) 



of vectors {v ^ m '}m=x, It is also easy to see that the matrix V with entries Vi :m = 
is symmetric. This is proved by noting that the generating function 



**HEA".!^4^, (2.8) 

is manifestly symmetric in x and y. From this symmetry and the orthogonality relation 
(p/7|) , we deduce the following completeness relation 



£ 4^ = • (2-9) 

m> 1 



The first statement above, that i>( m - ) be the m-th eigenvector of T(g,a), is proved 
analogously, by means of a contour integral over the unit circle: 



i>l ^ ' i i,j>l 



1 \m-l 



r j / dy gxy \{<l ~ y) 

V ~ Q f 2my (1 - gax) - y(ga + g*(l - a*)x) (1 - J )™ (2 ' 10) 

ga + g z (l — a z )x 1 — gax 

((1 - gax)q - {ga + g 2 (l- a 2 )x)) m ~ l 



1 — q 2 gx 



(l — gax — q(ga + g 2 (l — a 2 )x)) 



Note that only the pole y = yo(x) = (1 — gax)/(ga + g 2 (l — a 2 )x) has contributed to 
the contour integral, as we have taken now only the poles lying outside the unit disc, and 
yo(x) > 1 for |x| < 1, thanks to the inequalities ([O]). Using the equation (|2.5|) , we may 
write g 2 (l — a 2 ) = 1— ag(g+^), and notice that the numerator and denominator monomials 
in the last line of (|2.10| ) respectively read 

(l-gax)q- (ga + g 2 (l - a 2 )x) = (q-x)(l- — ), 

<1 (2.11) 

1 — gax — q(ga + g 2 (l — a 2 )x) = (1 — qx)(l — qga), 
so that we finally arrive at the eigenvalue equation 

T(g,a)v^ =A m v^\ (2.12) 

where 

(I — 9^\m-l 

A m = A m (g,a) = g - (r -±-^ r , (2.13) 

for m = 1, 2, 3, ... 

Let us examine our result in the pure Lorentzian gravity case a = 1. The quadratic 
equation ( |2.5| ) is solved as 

q(g,l) = gC(g 2 ), (2.14) 

for \g\ < i where 



cw = 1 - = y; W (2 . 15) 

2x 4^(n + l)!n! ' v ; 



n>0 

i2 1 



is the generating function of the Catalan numbers, satisfying xC (x) = C(x) — 1. The 
eigenvalues can then be simplified to read 

A m = [gCig 2 )] 2 ™- 1 = q 2m -\ (2.16) 



Note that as \q\ < 1, we have Ai > A2 > ... > 0, and all the eigenvalues have series 
expansions in powers of g with positive integer coefficients. 
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2.3. Commuting transfer matrices 

We now make the crucial observation that the eigenvectors v( m \g, a) of the previous 
section only depend on q, which is itself a certain function of g and a. Hence there exists an 
infinite family of matrices T(g,a) sharing the same eigenvectors, namely those for which 
the values of g, a lead to the same value of q. But matrices which can be diagonalized 
simultaneously form a commuting set. More precisely, the following statement holds: 

[T(g,a),T(g ,a ) = iff = — =q+~, (2.17) 

ga g'a! q 

for arbitrary q (assumed to be real and such that \q\ < 1). Moreover, the common eigenvec- 
tors are given by ( p?6|) . Again, Eq. (|2.17|) can be proved directly by use of contour integrals 



involving the product of two generating functions ( |2.3| ) . 

To better understand the mechanism of this commutation, let us express the transfer 
matrix T in terms of its orthonormal eigenvectors and of its eigenvalues, namely 



0{x,y\g : a) = ^ F m (x\q)k m F m (y\q) 
m>l 

= Y^EFmix^X^FUylq) (2 . 18) 



m>l 



m>l 



where we have identified 



I — 22. 

A=- q -, (2.19) 

1 - qga 



by use of ( |2.13| ). For any fixed q with \q\ < 1, we may use the (spectral) parameter A to 
characterize the different commuting matrices ( |2.17| ). Denoting 



T q (X)=T{g,a) for 1 ^ 2 ( 1 a ') = q + 1 and ga = J_A, ( 2 .20) 

ga q ± - qX 

and 

V q W=vW(g,a), (2.21) 
independently of A, we have from ( |2.18| ) 

T q (X) = VX X m ~ 1 V( rn \V q ( rn * > ) t , (2.22) 

m>l 
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where we have used the orthogonal projector = V q (Vq^Y onto the m-th eigenspace 
of T q (\). Thanks to the orthogonality relation for the eigenvectors, which translates into 
p( m ) ' p( m ) = S m>m i Pq m \ we easily get that 

T q (X)T q (X')=T q (XX'). (2.23) 

This relation trivially implies the commutation of transfer matrices fl2.17|) . But this is 
a much stronger constraint. Note that all real values of A and q such that \q\ < 1 are 
allowed. Note also that both the commutation (|2.17| ) and the multiplicativity property 
of spectral parameters ( |2.23| ) hold only with the help of the curvature weight a. The 
"solvability" of the case without a curvature term (allowing for the solution of |p.2| ) is 
simply a consequence of these more general properties. Finally, let us mention that the 
transfer matrix of the regular lattice is a member of the family indexed by the parameter 
q for any q, as it can be obtained by taking the limit o — > 0, g — > 1 while keeping the ratio 
(1 — g 2 (l — a 2 )) I (ga) = q + (1/q) fixed. In this limiting case, A = 1 and the transfer matrix 
T q (l) is nothing but the identity matrix. 

To conclude this section, let us rewrite the generating function ( |2.3|) in terms of q and 
A only 

e q (x,y\X)= 5>V[T 9 (A)] W 

(2.24) 



VXxy(l 



(1 - qx)(l - qy) - X(q - x)(q - y)' 

2.4- Squares and Triangles 

The curvature model introduced above may be reinterpreted as a discrete model 
for Lorentzian gravity, where the world-sheet is generated by arbitrary tessellations with 
squares and triangles that respect chronology. This is easily seen by performing the fol- 
lowing transformation on the configurations of the previous model. 

Let us expand again in powers of x and y the generating function 



Kx, y\g, a) = 9Xy (2.25) 

1 — ga(x + y) — g z (l — a z )xy 

and interpret each contribution to x l yi as a configuration of i lower and j upper half-edges. 
By expanding the denominator of ( |2.25| ), we get an arbitrary sequence of terms gax, gay 
or g 2 (l — a 2 )xy. Whenever a term gax (resp. gay) is picked, we interpret it as an isolated 
lower (resp. upper) half edge, which comes with a weight ga. Whenever a term g 2 (l — a 2 )xy 
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Fig. 3: The Boltzmann weights of the square and triangle formulation of the 
curvature model for discrete 2D Lorentzian gravity. We have also represented 
a typical world-sheet configuration in this new interpretation, together with 
its dual, made of squares and triangles. 

is picked, we now have a pair of a lower and an upper half-edge which we can regroup so 
as to form a crossing of the time line (see Fig. 3). Such crossings come with a weight factor 
g 2 (l — a 2 ). Finally, the numerator gxy in ( |2.25| ) is there to ensure staircase boundary 
conditions for these new configurations. It corresponds to adding a lower half-edge at the 
far left and an upper one at the far right, each of which comes with a factor yfg only (and 
no a). 

We thus get another representation for the configurations contributing to the transfer 
matrix element Tij(g,a), with the Boltzmann weights defined in Fig. 3. Dually, these 
correspond to supplementing the triangulations we have considered so far, by squares with 
two time-like and two space-like edges (see the example of a world-sheet configuration 
and its dual depicted in Fig. 3). All the weights have now been translated into different 
fugacities for the triangles and for the squares. Each triangle comes with a weight ga, 
while each square receives a weight g 2 (l — a 2 ). Note that this latter weight is positive 
or negative according to whether a is larger or smaller than 1. Finally, the case of the 
boundary triangles is special since they receive a weight y/g only. Let us also mention that, 
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in the limit a — > corresponding to a regular lattice, only squares and boundary triangles 
survive, in which case it is clear that the only possible arrangement for the tessellation is 
regular and the transfer matrix is the identity matrix as it should be. 
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Fig. 4: The connection between square-triangle tessellations and simple 
triangulations is done by splitting each square (i.e. each crossing) into two 
triangles (i.e into two half-edges). To recover the proper curvature weight, 
we transfer the a factors as shown by the arrows. We then add up weights 
corresponding to the same final triangulation. 

To make the contact with the curvature model, the idea is to split each square of 
the tessellation into a pair of triangles. To recover the proper curvature weights, we 
chose to split each square into a down triangle followed by an up one. Now the same 
triangulation is obtained from several square-triangle tessellations and the matrix element 
Tij(g, a) is obtained by summing over all these configurations. In the dual language, our 
splitting corresponds to separating each crossing into a lower and an upper half, such that 
the upper- half sits at the left of the lower one (see Fig. 4). As far as the g factors are 
concerned, since each crossing receives a factor g 2 , which after splitting, gives a weight g 
per half edge, it is clear that these factors are well taken into account. The yfg factors 
at the boundary are also correct. The a factors are more subtle. Before we collect all 
these factors, we remark that we can transfer the weight a of each lower half-edge, which 
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is not at the left boundary, to the horizontal link sitting just above it and to its left (see 
Fig. 4). Similarly, we transfer the factor a attached to each upper half-edge, which is not 
at the right boundary, to the horizontal link sitting just below to its right. Then all the 
horizontal links which separate two lower half-edges (resp. two upper half-edges) receive a 
factor a from their right lower half-edge (resp. their left upper half-edge) as they should in 
the curvature language, while all the horizontal links with a lower half-edge on the left and 
an upper half-edge on the right receive a factor 1 as they should. Finally the horizontal 
links with an upper half-edge on the left and a lower half-edge on the right receive a 
weight a 2 from the two half-edges, together with a contribution (1 — a 2 ) coming from the 
configuration where these two half-edges were connected to form a crossing. Adding these 
weights restores a factor 1, as required. 

2.5. Dimers 

From our experience with two-dimensional Euclidean quantum gravity we know that 
the simplest way of introducing matter in two-dimensional universes is to allow for dimers. 
Dimers are decorated links which have a certain weight associated with them and each 
triangle in a triangulation can carry at most one dimer. It is well known that by fine- 
tuning the weight of the dimers in two-dimensional Euclidean quantum gravity to some 
specific negative value one can reach a critical point where the universality class of the 
underlying geometrical system changes from that of pure gravity to that characteristic 
of gravity coupled to non-unitary matter with central charge c = —22/5. It is therefore 
interesting to note that our model of triangles and squares can be viewed as a simple dimer 
model. This interpretation simply comes about if one views each square as consisting of 
two triangles sharing a dimer which lies along either one of the diagonals of the square. 
In this picture the triangle building blocks do not have any dimers associated with them. 
The weight of a square, g 2 (l — a 2 ) is hence to be understood as 2 • {go) 2 ((1 — a 2 )/ (2a 2 )) 
where ga is the weight of a triangle, (1 — a 2 ) / (2a 2 ) is the weight of a dimer and the factor 2 
in front takes into account that there are two ways of placing the dimer inside the square. 
From this decomposition it follows that if new critical behaviour occurs it should happen 
at a value of a for which a > 1. However, it is important to note that our model only 
includes a subset of configurations of the full dimer model since we do not have any dimers 
on the space-like links. Thus, if we do not see any new critical behaviour for the model in 
question it does not necessarily imply that such behaviour does not occur for the full dimer 
model. In this connection, let us mention that there have been studies of two-dimensional 
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Euclidean quantum gravity coupled to dimers in a spirit similar to this one. Namely, one 
has calculated loop-loop correlation functions as a function of distance, not leaving out 
any matter configurations but modifying the concept of distance so as to avoid dimers on 
the entrance and exit loop. In that approach a continuum behaviour different from that 
of pure gravity was seen . 



2.6. Solution and continuum limit 

The property (2.23) yields immediately the generating function of the T-th power of 



the transfer matrix, describing a world-sheet of time-size T: 

6 q , T (x,y\\) = ^2 xl y Jz i,j( T \9,a) 
*,i>i 

= (2.26) 

»ii>i 

xy(l — q 2 )X^ 
(1 - qx)(l - qy) - X T (q - x)(q - y) ' 

for all T = 0,1,2,.... 

As noted in |TB|], one can derive a continuum limit of the expression ( |2.26| ), that 



corresponds in our notations to the limit q — > 1 (ga — > a/(a + 1)). The critical values 
x = y = 1 (corresponding to the limiting convergence radii of the series F m (x\q) and 
F m (y\q)) must be approached simultaneously. We see then from ( |2.26| ) that A must also 
tend to 1 for the expression to remain finite. But from ( |2.19| ), a natural way to realize the 
latter limit is to simply keep a constant while q — > 1 (and g — > l/(a + 1)). 

Let us therefore perform the following scaling transformations for some small param- 
eter a 

T=-, x = l-aX, y = l-aY, q = e - a " /X , (2.27) 
a 

where A is the renormalized cosmo logical constant, r the continuous time variable, and a 
an infinitesimal parameter with the dimension of a length. This corresponds in turn to the 
scaling behaviour 

A = 1 -2aav / A + C(a 2 ), g = — ^(1 - ^« 2 A) + C(a 4 ), (2.28) 

for fixed a. Applying the transformations Q2.271) to ( |2.26|) , we get the rescaled two-loop 
generating function 

6(A, Y\t, VX, a) = lim a9 qT (x, y\X) 

2VAe~™^ (2-29) 



(X + y/X)(Y + VX) - e- 2Ta ^(X - VX)(Y - VX) 
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The actual two-loop correlator 

G{L u L 2 \T,yfK,a) = lim i-Z hM (T\g, a), (2.30) 



with = ali, i = 1, 2 is obtained by performing the inverse Laplace transformation of 
( |2.29|) , leading to (we have set <fi = e~ aT ^ for simplicity) 

/too 
dAWe LlX+L2Y 6(X, F|r, v 7 !, a) 
-ioo 

joo 2v / A(/ , e v/ ^^ LlU+L2,; ^ 
. loo dudV (u + l)(v + l)-^(u-l)(v-l) 

rioo e VA^Li^-L 2u(l _^, ) + (l+ ^ 2) 



/zoo 
«(l-0 2 ) + (l + 

Va( L x u 

j_ju — _ ^ - - — (2 . 31) 



VX( L 1U +L 2 ^ u(1 _^ ) 1 +(1+ ^2 ) 



2VA0e-^ L2 ^^ 



fc>0 



k\ J- loo Ml - </> 2 ) + (1 + 



l-</> 2 ^ (A;!) 2 V(l-</> 2 ) 2 7 



i-<#> 2 



where we have first made the change of variables X = u\/X and Y = vy/X, and then 
used the Cauchy formula to express the integral as a sum over residues (only one pole in 
v contributes, located at v(u) = — (it(l + (j) 2 ) + (1 — (j) 2 ))/(u(l — <p 2 ) + (1 + 4> 2 )), whereas 
only the multiple poles at u = (1 + (p 2 )/(4> 2 — 1) occur). In terms of the rescaled variables, 
the two-loop correlator reads finally 

G(L 1 ,L 2 \r, \/X,a) 

( \ (2-32) 



-V / A(Li+L 2 )cotanh(arV / A)j o ( ^\/^-L 1 L 2 \ 

\sinh(arV / A) / 



sinh(arV / A) \sinh(aT\/A) 

where Jo is the modified Bessel function Io(x) = Y^k>o( x /^) / i^) 2 ■ 

A few remarks are in order. First note that the dependence on the curvature parameter 
a is quite simple, and that the physics of the model is not affected by it. So we must draw 
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the rather negative conclusion that the introduction of curvature in the model cannot 
change its scaling behaviour. It has just shifted the critical value of g (from 1/2 for the pure 
gravity case to l/(a + 1) in the model with curvature). Moreover the explicit dependence 
on a can be entirely absorbed in redefinitions of Li,L 2 ,A and G, namely L\ — > aL\, 
L 2 — > aL 2 , A — > A/a 2 and G — > aG. This "triviality" of the curvature dependence will be 
explained in Sect. 3 below, from a very different perspective. 

As a final consistency check, let us compute the a — > limit of the two-loop correlator 
( p.32| ). Using the asymptotics of the Bessel function Iq(x) ~ e x I \J1txx for large x, we get 



G(L u L 2 \r, VA,a-^0) 



V AnTa\fLiL 2 (2.33) 
' 6(y/L^- yfT 2 ) = 8{L 1 -L 2 ), 



2JL 



where we have used the limit lim s ^ + e x l s j^pKS = 8(x). We recover therefore the 
expected flat world-sheet result (also stated as Zi lt £ 2 (T\l,0) = $i lt £ 2 )- 

2. 7. Other Boundary conditions 



In [13], where boundary conditions were chosen to be periodic, the result for the two- 



loop correlator in the scaling limit is very similar to ( 2.31 ) except that a = 1 (no curvature) 
and more importantly that the modified Bessel function Iq is replaced by I\. We have a 
very simple explanation for this fact, that involves computing the p-seam loop correlator 
which we will define now. Up to now, we only considered triangulations with "staircase" 
boundary conditions. Let us denote these triangulations as being of type (I). Starting 
from such a type (I) triangulation in its dual picture, a seemingly easy way to get periodic 
boundary conditions would be simply to identify the left and right staircase boundaries of 
the triangulation to construct a cylinder marked by a seam, remnant of the staircase. This 
construction is however problematic. Indeed, on each horizontal step of the left staircase 
boundary, we can have an arbitrary number of incoming lower half edges. Similarly, on any 
horizontal step of the right staircase boundary, we have an arbitrary number of incoming 
upper half-edges. When connecting these two staircases, there is no well defined natural 
prescription for deciding how to place these lower and upper half edges with respect to 
one-another. 

To overcome this problem , we need to introduce a new type of triangulations, which 
we will call type (II) , and which can be glued without ambiguity to a type (I) triangulation. 
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(I) 



(II) 




Fig. 5: An example of triangulation of type (I) and a triangulation of type 
(II). The boundary chaplets serve as seams to connect type (I) to type (II) 
triangulations. We give an example of cylinder made of two glued triangula- 
tions. Any marked point (black dot) on the lower (resp. upper) loop defines 
a unique chaplet of triangles pointing up (resp. pointing down). 

Let us first see how the existence of a "staircase" in the dual representation translates into 
the original (i.e. non-dual) representation made of triangles. Each leftmost lower half- 
edge of the dual picture becomes a triangle pointing up. Each such leftmost triangle is 
attached by its lower left vertex to the upper vertex of the leftmost triangle in the layer 
just below, so as to form a chaplet (see Fig. 5). Similarly, the right staircase translates into 
a chaplet of triangles pointing down, each triangle being attached by its upper right vertex 
to the lower vertex of the triangle in the layer just above. The "space" between these two 
chaplets is filled with strips of triangles of arbitrary length, including strips of length zero 
corresponding to the case where the two chaplets would be in contact. The two chaplets 
will serve as seams in our construction. In order to glue our type (I) triangulation to a 
type (II) triangulation, this type (II) triangulation must itself have boundaries made of 
chaplets but now with the reverse convention, i.e. with a chaplet of triangles pointing up 
as its right boundary and a chaplet of triangles pointing down as its left boundary (see 
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Fig. 5). Again, the space between these chaplets is to be filled with strips of arbitrary, 
possibly zero, length. It is clear that we can glue a type (I) triangulation to a type (II) 
triangulation on any side by simply superimposing the left (resp. right) chaplet of the 
former to the right (resp. left) chaplet of the latter. We can thus construct a p-seam loop 
correlator by gluing p triangulations of alternating type (I) and (II), and, assuming that 
p is even, gluing the last triangulation to the first one so as to form a cylinder (see Fig. 5 
for p = 2). If the number p of triangulations is odd, we can build an open object with 
p — 1 seams and either type (I) or type (II) boundary conditions, depending on whether 
to two extremal triangulations are of type (I) or of type (II). This construction gives a 
nice a posteriori explanation for our choice of weight yfg per boundary half-edge: this is 
because the edges (or the triangles) have to be identified by pairs to form the seam, hence 
the weight is right to get a factor g per half-edge of the seam. 




(c) (d) 

Fig. 6: The one-to-one correspondence between type (I) and type (II) trian- 
gulations. Starting from a type (I) triangulation (a), we move each slice one 
step to the left with respect to the slice just below to obtain (b). We then 
send the lower left and upper right vertices of the boundary triangles thus 
released in each slice to respectively upper left and lower right position to 
obtain the type (II) triangulation (c), or equivalently (d). The passage from 
(a) to (d) is clearly invert ible. 

We have computed in Sect. 2.1 the transfer matrix for type (I) triangulations, given 
by (|2.2|) . A remarkable result is that, although type (II) triangulations are different from 
type (I) triangulations, the transfer matrix for type (II) triangulations is identical to that 
of type (I), i.e. 

T$\g, a) = T$(g, a) = T^g, a), (2.34) 

with Tij(g, a) given by fl2.2p . To understand this property, we remark that starting from 
a type (I) triangulation, we can deform it into a type (II) triangulation by simply letting 
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each time slice move one step to the left with respect to the slice just below, and by sending 
the lower left and upper right vertices of the boundary triangles so released in each slice, to 
respectively upper left and lower right positions (see Fig. 6). This transformation is clearly 
invertible, hence the announced result ( |2.34|) . 

Thanks to the above construction, we may define and compute the p-seam two-loop 
correlator by juxtaposing side by side p copies of alternating (I) and (II) triangulations, 
with a total lower (resp. upper) loop length of l\ (resp. £2), in which we simply have 
to identify the consecutive chaplet boundaries (again, the boundary weights produce the 
desired weight g per half-edge of seam) . Assuming p even and gluing the last triangulation 
to the first one, this leads to the p-seam two-loop correlator 

Z^ iia (T\g, a) = [(T(g, a) ® T(g, a) ® ... ® T(g, a))\^ (2.35) 

where there are p factors in the tensor product, and the matrix element (£±,£2) actually 
corresponds to the sum over all p-tuples of pairs of indices (raj, rii) for each T(g, a), with 
Y^ii n i ~~ 1/2) = £\ and Yli( m i — 1/2) = ^2- In terms of generating functions, the p-seam 
correlator (|2.35|) is simply generated by: 

= (^^Y ', peven. (2.36) 

v v x y j 

For odd p, gluing p triangulations of alternating type (I) and (II) leads to an open object 
of either type (I) or type (II) with the two-loop correlator now generated by: 

^(^y|A) = v^f ^ T( !l5 |A) y, podd. (2.37) 

Let us now perform the scaling transformations (|2.27[) on ( |2.36| ) or ( |2.37|) .We obtain 
the p-seam scaling function (note that it has required a multiplication by a p to produce a 
finite limit) 

6 (p) (X, Y\t, v 7 !, a) = G(X, Y\t : VX : a) p . (2.38) 
We get the corresponding two-loop correlator by applying the inverse Laplace transform, 
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and performing integrations analogous to those of (|2.31 ): 

/ioo 
dXdYe LlX+L2Y Q(X, Y\t, VX, a) p 
-ioo 

= A r dudve^ u+L ^ ( -= ^ Y 

J-ioo \VA((u + l)(v + 1) - <P 2 (u - l)(v - I)) J 



(p-1)! (VX(u(l - <\> 2 ) + (1 + <P 2 )))p dvP- 1 6 



VX(2(/>L 2 )p f io ° , e 
du 



VAI Liu-L 2u(i _^ ) + (i + ^ 2) 



(p-l)ILa J_ loo («(l-02) + (l + 02))p 



(p - 1)!L 2 ^ fc! Vl - </> 2 7 («(1 - </> 2 ) + (1 + </> 2 )) p+fc 

= 2</>y/A( v /L^)P- 1 _VA (Ll+L2) i±^ v (jrgVALiLaj 
(l-^)(p-l)! 6 ^ fc!(* + p-l)! 

(2.39) 

for p > 1. In terms of the rescaled variables, the p-seam two- loop correlator reads 

VXiVL^r- 1 „-v^a 1+i9 kotanhfarVAU ( I^KL^ \ (2-40) 



sinh(ar\/A)(p — 1) 



-e 
I 



-V / A(Li+L 2 )cotanh(arVA)j ( 2VAL1L2 \ 

P \sinh(arv / A) / 



where I m denotes the m-th modified Bessel function, defined by the series I m (x) = 

When p = 1 this reduces to (|2.32|) , as it should. More interestingly, when p = 2, in 
the pure gravity case a = 1, the expression reduces, up to a factor of L2, to the result 
for the two-loop correlator in periodic boundary conditions with one marked point on the 
lower loop [[12], i.e. equivalently, the two-loop correlator in periodic boundary conditions 
with one marked point on the lower loop and one marked point on the upper loop (note 
that marking a point on the external loops simply amounts to multiplying by the length 
of this loop in the correlator). 

On can easily explain this coincidence by noticing that, starting from a cylindric 
world-sheet with a marked point, say, on the lower loop, this marked point defines a unique 
chaplet of up triangles crossing the cylinder from the lower to the upper loop. Indeed this 
point is the lower left vertex of a unique up triangle in the first slice, whose upper vertex 
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is itself the lower left vertex of a unique triangle in the second slice, and so on (see Fig. 5). 
Similarly, any marked point on the upper loop defines a unique chaplet of triangles pointing 
down. Thus marking a point on the upper loop and one on the lower loop amounts to 
marking two chaplets, one made of up triangles and one made of down triangles. Such 
chaplets cannot intersect and thus divide the cylinder into two triangulations, one of type 
(I) and one of type (II). This explains the connection between the two-seam two-loop 
correlator and the periodic loop correlator. Note that this nice property breaks down as 
soon as p > 2 because, for more than one point on each external loop, we cannot guaranty 
for arbitrary chosen marked points that the corresponding chaplets alternate along the 
cylinder between "up" and "down" chaplets, which is crucial to keep our interpretation as 
p-seam correlator. In particular two successive chaplets of the same (up or down) type can 
now merge into a single chaplet. 

Note finally that we have the following relation between our (2p + 2)-seam correlator 
and the amplitude A p (with p + 1 a type of winding number) obtained by Nakayama in a 
continuum calculation using proper time gauge JT1J 



3. Lorentzian Triangulations as Random Walks 

In this section, we will discuss the equivalence between Lorentzian gravity configura- 
tions and ordinary random walk (RW) configurations in one dimension. This connection 
will allow us to re-phrase the quantities computed above in the language of RW statis- 
tics. For instance we will see how the two-loop correlator of Lorentzian gravity relates 
to a well-known generating function for the large excursions of a Brownian motion. This 
equivalence will also explain why Lorentzian triangulations have fractal dimension 2 and 
why introducing a curvature weight (i.e. taking 1) cannot change the continuum large 
scale properties of Lorentzian gravity, as found in Sect. 2.4. Finally, the RW picture will 
also provide a very direct derivation of the crucial property (|2.23| ) for the product of two 
transfer matrices. 




(2.41) 
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3.1. Equivalence between Lorentzian triangulations and random walks 

The random triangulations describing Lorentzian gravity can be seen as deformations 
of a regular triangular lattice. More precisely, if we view the regular triangular lattice as 
made of regular time slices of alternating up and down triangles, then removing some of 
these triangles and gluing together the remaining ones along their time-like edges within 
each strip leads to a triangulation like that of Fig.l. Conversely, we can, starting from 
a random triangulation of this type, insert a number of additional triangles to make it 
regular. It is therefore possible to visualize a random triangulation as a regular one with 
two types of triangles: "real" triangles which are to be kept and form the triangulation 
and "virtual" triangles which have to be shrunk or removed. 

This splitting procedure into real and virtual triangles can be made according to a well 
defined procedure which we shall now describe, and which takes the form of a one-to-one 
correspondence between Lorentzian triangulations and directed random walks drawn on 
the regular triangular lattice. 

The equivalence is as follows: 

Consider a directed random walk drawn on the regular lattice, starting at "height" 
h = 0, making Ah = ±1 steps to the right and ending at height h = T + 1 after S steps 
(see Fig. 7). Since the walk is directed, the horizontal direction can be viewed as well 
as an effective time dimension for the random walk whose motion takes place in the one 
dimensional vertical /i-direction. This height variable h will however correspond to the 
time variable t of the triangulation, with each time slice t lying between heights h = t — 1 
and h = t. The random walk is moreover required to stay confined within the strip 
< h < T + 1. Now to each elementary ascending step of the RW between heights h and 
h + 1, we associate the pair of triangles of the regular lattice made of the up triangle lying 
immediately to its right in the time slice h + 1 and the down triangle lying just below 
in the slice h (see Fig. 7). These triangles will be the "real" triangles to be kept in the 
triangulation. Removing all the other triangles and gluing the real ones together along 
their time-like edges, we end up with a random Lorentzian triangulation of width T, and 
decorated by l\ down triangles in time slice t = and I2 up triangles in time slice t = T+ 1. 
These extremal time slices are then removed to recover the relevant triangulation. Note 
that from the above construction where we add triangles only to the right of the ascending 
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(b) (c) 

Fig. 7: Starting from a directed random walk (a) drawn on the regular 
triangular lattice, we associate to each ascending step of the walk the pair of 
triangles lying just below it (grey triangles) . We then eliminate all the other 
triangles and glue the grey triangles together along their time-like edges inside 
each strip, thus forming a random Lorentzian triangulation (b). The and 
T + l slices (light-grey triangles) are finally removed (c). 

steps, it follows that the resulting triangulation satisfies the staircase condition of Sect. 
2.10. 

We can thus associate to each directed RW of the type described above a Lorentzian 
triangulation with staircase boundary conditions. Let us now describe how, starting now 
from such a Lorentzian triangulation, we recover the associated RW. This reverse construc- 
tion can be performed in three steps as shown in Fig. 8. The first step is to go to the dual 
representation of the triangulation (Fig.8-b). In a second step, the dual world-sheet con- 
figuration can be transformed into a tree configuration (Fig.8-c) by attaching each vertical 
dual bond to the bond sitting in the strip just below and to its left. In a third step, this 
tree-like structure is transformed into a directed random walk (Fig.8-d) by simply following 

7 Note that the triangles can be viewed as well as being added on the left side of the descending 
slopes of a directed walk, now going from height h = T and reaching height h = — 1 after S steps 
to the left. 
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Fig. 8: Starting from a Lorentzian triangulation satisfying staircase bound- 
ary conditions (a) , we go to the dual configuration (b) and attach each vertical 
dual bond to the bond to its left in the preceding strip, as shown by the ar- 
rows, to get a tree structure (c). Following the contour of the trees from 
the lower left branch to the upper right one (black points), the sequence of 
ascents and descents can be translated into a directed random walk (d). 

the contour of the trees from the lower left branch to the upper right one and making a 
step Ah = +1 for each ascent along a branch and a step Ah = — 1 for each descent along 
a branch. It is easy to check that the resulting RW is precisely the expected one. 

Our main result here is thus a one-to-one correspondence between directed random 
walks and Lorentzian triangulations, which in turn implies a direct equivalence for the 
various quantities describing them. Some of these relations will be described in the next 
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section. Let us simply note here that we can deal with all types of boundary conditions. 
The staircase boundary condition described in this section corresponds to random walks 
starting at h = and ending at height h = T + 1. The case of walks starting and ending 
at h = would describe periodic boundary conditions with a marked point in the lower 
loop of the triangulation. 

3.2. Two-loop correlator as a generating function for large excursions 

Let us now re-interpret some of the natural statistical properties of random Lorentzian 
triangulations in the language of random walks. 

Since we associate a pair of triangles to each ascending step, the number N t of triangles 
is directly related to the length £ of the walk by: 

N t + (h + l 2 )=S + T + l. (3.1) 

Here N t counts all the triangles with 1 < t < T, ignoring the (l\ + l 2 ) triangles in the 
slices t = and t = T + 1 which are removed in the construction of the triangulation. We 
also used the fact that there are T + 1 more ascending steps than descending steps in the 
walk. Setting a = 1 for the time-being, the triangulation comes with a factor g Nt ~ T due 
to the particular treatment of the extremal left and right triangles in each of the T slices. 
The corresponding weight g s -(h+ l 2- 1 ) can De interpreted as a fugacity g per step for the 
random walk, together with a reflection factor (1/g) for each of the l\ contacts at h = 
or the l 2 contacts at h = T + 1, except for the last contact. Note that for the critical value 
g = 1/2, we can view these weights as a probability 1/2 for the walk to make a step up or 
down inside the strip and a probability 1 for the walk to make an up (resp. a down) step 
at the h = (resp. h = T + 1) boundary, i.e. we have an unbiased random walk between 
two reflecting walls. 

The sizes l\ and I2 of the bottom and top loops of the triangulation correspond pre- 
cisely to the number of contacts at h = and at h = T + 1 respectively. We can thus 
reinterpret the two-loop correlator Zj lj j 2 (T|y, a = 1) as the generating function 

Z luh (T\g,a=l) = 1 J2( 2 9) Sp T(h,h,S) (3.2) 
[A9) s>o 

where Pr(hj S) is the probability for an unbiased walker starting at h = and evolving 
between two reflecting walls at h = and h = T + 1 to reach the wall h = T + 1 in S steps 
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after exactly li contacts (including the first one) at the bottom wall and I2 (including the 
last one) at the top wall. 

More generally, if we re-instate the factor x (resp. y) per space like link in the bottom 
(resp. top) loop, the triangulation generating function 9t(x, y\g, a = 1) (given by Eq. ( |2.26p 
for a = 1, i.e. q = (1 — a/1 — 4g 2 )/(2g) and A = q 2 ) becomes the generating function for 
random walks between two walls with a weight g per step and a weight x/g (resp. y/g) per 
contact to the bottom (resp. the top) wall. In the RW language, we call each portion of 
walk between two successive contacts an "excursion" . Since the walk is a Markov process, 
each of these excursions is an independent random object. Excursions here are of four 
types: bottom-bottom and top-top excursions, which have the same generating function 
(up to the exchange of x and ?/), and bottom-top and top-bottom excursions which also 
have the same generating function. For instance the bottom-bottom excursion counts the 
walks starting and ending at h = and staying below h = T. The distribution probability 
for, say, the height and the length of bottom-bottom excursions is a basic quantity of 
the random walk statistics and enters many different physical problems. Its continuous, 
Brownian motion limiting scaling expression is well known, as are many other more involved 
excursion properties appearing in many different areas of physics. As an example, we would 
like to show here how the two-loop correlator of Lorentzian gravity directly relates to a well 
known probability distribution describing the diffusion of a particle in a one-dimensional 
random energy landscape. 

The diffusion of a particle in a random energy landscape can be modelled by the so- 
called Sinai model ||20|1 , where a particle jumps to one of its two neighbouring sites with 
a probability p oc exp(— AU) involving the energy difference AU with this neighbouring 
site. If we represent the energy landscape as a random walk, with the height h oc U 
corresponding to the energy, a bottom-bottom excursion describes an energy barrier of 
height less than T. Exact expressions for the large time properties of the diffusion have 
been obtained and it was recognized that a good way to recover these expressions is to 
assume that the particle at time t is localized precisely in the lowest energy minimum it 
can reach by passing all the energy barriers of height less than a certain maximal size Tit), 
with the relation T{t) oc ln(t). Let us for instance consider the example of Fig. 9 (b) where 
a particle, starting at a wall, jumps to the right by passing all barriers of height strictly 
less than T, looking for lower and lower minima to the right until it reaches a barrier of 
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T+l 




Fig. 9: (a) A random walk entering the definition of Z 7n+ i i i(T|(7, a), i.e. 
with m + 1 contacts (here m = 4) at h = and 1 contact at h = T + 
1. (b) The corresponding energy landscape obtained by erasing in (a) all 
the ascending steps following immediately a contact at h = 0. The deepest 
minimum reachable by a diffusing particle (black dot) passing energy barriers 
of height strictly less than T is at depth m (white circle) . 



height larger than or equal to T. How far and how deep does the particle go? From ( |3.2| 
we know that the quantity 



Zi 1=m+ i^ 2= i(T\g, a = 1) 



1 



771 + 



T ^(2^) s P T (m + l,l,^) 



(3.3) 



S>0 



relates to the probability PT{m + 1, 1, S) that the walk reaches h = T + 1 for the first time 
after S steps and passes m times at h = 0. If we deform the walk by suppressing the first 
ascending step after each contact point at h = 0, the quantity Prim + 1,1, S) becomes 
the probability for the diffusing particle to have reached a depth m by passing barriers of 
length strictly less than T and being blocked after S — m — 1 steps by a barrier of height 
larger or equal to T. Going to the continuous limit and using ( |2.32| ), we obtain that the 
two-loop correlator 



G(L,0\t = 1, VX,a = 1) 



should be equal to the Laplace transform 



sinh(v / A) 



-\/ALcotanh(\/A) 



(3.4) 



P(L, A) = / dxP(L,x)e 
'o 



-Ax 



(3.5) 
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of the probability density P(L : x) for the particle to reach the deepest minimum at depth 
L and being blocked at a distance x from the wall, in the scaling limit of large T, m, and S 



Pi(L,xi) that the particle reaches the deepest minimum at depth L and distance x\ from 
the wall, and the probability ^2(^2) that the last excursion to the blocking point is of 
length x 2 - The expression (|3.4|) is then a direct consequence of the following expressions 



which are well known in the context of the Sinai model |21j] . 

The RW equivalence also provides an explanation for the fractal dimension 2 of 
Lorentzian triangulations. It is well-known that the Brownian motion has fractal dimen- 
sion 2, meaning that the time (represented here by the length S of the walk) needed to 
travel over a distance T (here the height h ~ T) scales as T 2 . Since S corresponds in our 
equivalence to the number N t of triangles, i.e to the area A of the triangulation, we obtain 
that A ~ T 2 and thus get a fractal dimension 2 for Lorentzian triangulations. Note that 
the fact that l\ (or I2) scales as T is also clear in the Sinai model language (with l\ ~ m) 
where it means that the depth reached by passing maximal barriers of order T is itself of 
order T. 

3.3. Curvature and universality 

Let us now re-instate a curvature weight a 7^ 1 in the problem, and let us see how 
it translates into the random walk language. As for the equivalence with tessellations of 
triangles and squares of Sect. 2.4, we can regroup the up and down triangles in each strip 
as follows. One first re-groups all pairs of triangles made of an up triangle immediately 
followed to its right by a down triangle and assign a weight g 2 to these pairs. The remaining 
up triangles necessarily have an up triangle to their right. We thus assign them a weight ga 
accounting for the fugacity g per triangle and the curvature weight a for their interaction 
with their right neighbour. Similarly, each remaining down triangle necessarily has a 
down triangle to its left and again receives a weight ga for the same reasons. With this 

8 Strictly speaking, the particle is blocked after S — m — 1 steps, not S, but this correction 
scales less rapidly than T 2 and is thus negligible. 



with m/T = L and S/T 2 



x fixedtl. This probability is the convolution of the probability 




— VAL cot anh ( ) 



(3.6) 
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construction, all the weights have been collected properly and we simply need an extra 
factor 1/g to correct for the special boundary weights. Otherwise stated, we can assign 
a weight ga for each triangle and correct with a factor 1/a 2 for each pair made of an 
up triangle followed by a down one and a global 1/g factor. This way of correcting is 
alternative to that of Sect. 2.4 consisting in adding squares in the game. Denoting by P 
the total number of up-down pairs in the whole triangulation, we get a weight: 

1 ^ (3.7) 



9 (a 2 . 

Now in the RW picture, it is easy to see that up triangles lie immediately to the right of 
an ascending step, and that they are followed by a down triangle if and only if the next 
step is also ascending. In other words, pairs of up-down triangles correspond exactly to 

pairs made of two successive ascending steps. Denoting by n ++ (resp. n_| , n |_, n__) the 

pairs of successive ascending-ascending (resp. ascending-descending, descending-ascending, 
descending-descending) steps, we have the relations 

n+- = n_ + , 
P = n ++ , 

(3.8) 

T = n++ — n__, 



S = n++ + n_| + n |_ + n__ . 



The first relation simply states that there are as many summits as valleys in the RW 
landscape. Thanks to these relations and to Eq. (|Q|), we can rewrite ( |3.7|) as 



S— (h+h — 1) a n+++n-\ — +n — |-+n__+T— (h+h — 1) g 1 

a-+++---+ T = 9 ( ga ) ll+h - ian+ ~ ' (3 ' 9) 



The introduction of a ^ 1 has thus two effects in the RW weight: (i) it changes the 
interaction of contact with the walls from (1/g) to (1/ ga) and (ii) introduces a factor o for 
each change of slope in the walk. Note that the two effects cancel exactly at the walls so 
that the net contribution is simply a factor a for each change of slope inside the strip. Such 
a weight is nothing but an extrinsic curvature energy E oc — ln(o) for the one dimensional 
random walk. It is thus remarkable that the intrinsic curvature weight of two-dimensional 
Lorentzian gravity translates into an extrinsic one-dimensional curvature weight for the 
equivalent random walk problem. Note that according to whether a is less or larger than 
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one, one has a positive (favoring straight lines) or negative (favoring U-turns) curvature 
elastic constant. 

It is well known that extrinsic curvature is an irrelevant perturbation for the large 
scale properties of the random walk. At most it introduces a finite correlation length below 
which the random walk is rigid (positive curvature) or crumpled (negative curvature) but at 
scales larger than this correlation length, the curvature plays no role. Up to an appropriate 
rescaling, one thus recovers asymptotically the same continuum limit as without curvature 
energy. This is precisely what we observed in ( |2.32|) . 

3.4- Product of two transfer matrices 

From the RW picture of Lorentzian triangulations, we can easily re-derive the relation 
( 2.23Q for the product of two transfer matrices. Indeed, a two-matrix product element 
like [T(g, a)T(g', a')]ij corresponds to a sum with appropriate weights over all Lorentzian 
triangulations made of exactly two time slices, with a fixed number i of up triangles in 
their lower slice, a fixed number j of down triangles in their upper slice, and with staircase 
boundary conditions. In the RW representation, we thus have to consider all the random 
walks in a strip of height three < h < 3, starting at h = 0, ending at h = 3 and having 
exactly i contacts at h = and j contacts at h = 3. We thus have to sum over all these 
RW with a fixed number of contacts at the top and at the bottom of the strip < h < 3. 

Now we can decompose any such random walk into i + j — 1 "blocks" corresponding 
to the i + j — 1 portions of the walk between two successive contacts (see Fig. 10). The 
blocks are of four types, numbered (1) to (4), according to whether these contacts are 
at the bottom or the top of the strip, namely (see Fig. 10): (1) bottom-bottom, (2) top- 
top, (3) bottom-top, and (4) top-bottom. These blocks differ only by the (up or down) 
nature of their first and last steps while the intermediate steps form a saw-tooth sequence 
Ah = a, —a, a, —a, a, — a, ■ • ■, with a = 1 or — 1, of arbitrary (possibly zero) length k. The 
desired sum over the random walks can be achieved by first choosing one of the ( l ~^_^ 2 ) 
allowed sequences for the bottom and top contacts, then summing over all the random 
walks having this particular sequence of contacts, and finally summing over all sequences. 
Alternatively, we can attach directly a global weight individually to each of the four block 
types by summing over all possible intermediate configurations, i.e. over all values of the 
lengths k. 

If we are interested in the matrix-product T(g,a)T(g' , a'), we must attach a weight 
ga to each ascending step in the slice < h < 1, a weight g'a' to each ascending step in 
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(4) 




Fig. 10: The decomposition of a walk of height T + 1 = 3 into blocks. Each 
block is a portion of the walk between two successive contacts at the walls 
h = or h = 3. The blocs are of four types according to the four possibilities 
(1) bottom-bottom, (2) top-top, (3) bottom-top and (4) top-bottom for the 
two contacts at their extremities. 

the slice 2 < h < 3 and a weight (gg'aa') to each ascending step in the slice 1 < h < 2. We 
then must correct this weight by a curvature factor 1/a 2 for each sequence h = 0, 1, 2 and 
a factor 1/a' 2 for each sequence h = 1, 2, 3. Gathering these factors and summing over k, 
we obtain the block weights: 



w 



(i) 



(2) 



ga-^, . h. ga + g 2 (l — a 2 )g'a' 

9 a +h 22(99 'aa') k = y - f , , , 

or f— f 1 — qq'aa' 

a' 2 )ga 



k=i 



, , , g'a' 

9 a + — 

a' z 



w 



w 



(3) 



k=i 



99 

g'a' + g' 2 (l 



1 - sr^r'aa' 



9 2 9' 2 



(3.10) 



a^a 



fc=i 



gg'aa' 



(4) 



oo 

= V [gg'aa'f = . 

' 1 - gg'aa' 
fc=0 yy 



Note that all the weights have been gathered inside the blocks and that consecutive blocks 
do not interact. In order to have the commutation T{g, a)T(g', a') = T(g', a')T(g, a), it is 
sufficient that the four weights above are invariant under the change (g, a) <-> a'). This 



32 



is readily the case for the weights w% and while w% and are exchanged. The 
commutation condition reads therefore: 

defining for the two couples (#, a) and (<?', a') a common value of q as in (|2.5|). 

To go further and show (|2.23| ), we must compare the product T(g,a)T(g',a') to a 
single matrix say T(g", a"). Let us note that the same block decomposition as before can 
be achieved for random walks in a strip of height two, describing one-slice triangulations 
contributing to the elements of a single matrix T(g", a"). In this case, the contact sequence 
fully specifies the random walk and the four block types correspond the the four possible 
two step sequences: (1) Ah = +1,-1, (2) Ah = -1,+1, (3) Ah = +1, +1 and (A) Ah = 
— 1,-1. This leads immediately to the following block weights: 

wP = wf ] = g"a", 

™i 3) = g" 2 , (3-12) 

«4 4) = i. 

Let us now consider the same sequence of blocks with blocks of type (i), in both 
the product T(g, a)T(g' , a') and T(g",a"). On one hand, the weight contribution to 
T(g,a)T(g',a') is: 

-L («?>)" («£>)" («£>)" (4 4) ) , (3.13) 

where the pre-factor 1/gg' comes from the special weight of the boundary triangles. On 
the other hand, the contribution to T(g",a") is simply: 
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-(wV) U 3) ) . (3.14) 



Using the commutation relation ( |3.11|) and the fact that n*- 3 ** = + 1, i.e. the number 



of "ascending" blocks is one unit more than the number of "descending" blocks for the 
random walks that we consider, we can rewrite ( |3.13|) as: 

-1 . . „(l)_l_„(2) , s „(3) 

| / (l)\ n +n f (3) (4)\ n , 1R v 

wy) (wywy) . (3.15) 



/ ( 4 ) 

99 w 2 
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We get therefore obtain the product relation: 

T(g, a)T(g', a') = -^--T{g" , a"), (3.16) 
99 ' w 2 

provided that (i) the commutation relation ( p.llp is satisfied and (ii) we have the relation: 

(!) (!) // 99' 

21 y 1 - gg aa 



ga + go! - gg ad I q + 



(3) (4) (3) „ \ Q 

Wo Wr, = w\ a = 7 * - 



1\ (3.17) 



99 

with q as in ( |3.11| ) above. It is straightforward to check that the above values of g" and 
a" obey the relation: 

l-a" 2 (l-a" 2 ) 1 

9 , „ ] =9+-, (3-18) 
g a q 

and thus the three couples (g,a),(g',a') and (g",a") correspond to the same value of q. 
If we introduce the parameter A as in ( |2.19|) for (g,a) and A' for (g',a'), the value of 
A" for (g",a") is found to be A" = AA', while the pre-factor in ( |3.16| ) disappears since 
g" /gg'w^ = 1- We have therefore proven the relation ( |2.23j ) in the language of random 
walks. One should note that the power of this method comes from the fact that we used 
only four basic objects, the four blocks, each of which corresponds to a particularly simple 
family of walks, while the arbitrary sequence of the blocks plays no role. We expect that we 
can take advantage of such a block decomposition to obtain commutation relations for the 
transfer matrices of more involved problems of Lorentzian gravity, in particular including 
matter degrees of freedom. 



4. Other Models 

4-1. General construction: symmetric case 

Let us reconsider the expression for the curvature transfer matrix T q (X) in terms 

of the generating functions F m of its eigenvectors. If we now start from some arbitrary 
orthonormal family v^ m \ m = 1,2,... of vectors generating some space E, and a given 
non-zero function /(A), we may easily construct a commuting set of symmetric transfer 
matrices T(A) with entries 

^•(A) = /(A)5>i m) A— (4.1) 

m>l 
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that satisfy the multiplicativity property 

T(A)T(A') = ^^P^(AA'). (4.2) 

We also remark that, given any set of orthonormal vectors, generated say by F m (x), 
then the other set generated by G m (x) = F m (x)(p(x) is also orthonormal, provided we take 

k 

q — x 



<p(x) = T\(p(x\ qj ) k i , <p{x\q) = f , (4.3) 

xx 1 — qx 

3=1 



in which qi, ...,qk are arbitrary complex numbers inside the unit disc and kj some given pos- 
itive integers. Indeed, the functions <p(x\q) satisfy the inversion relation <p(x\q)<p(l/x\q) = 
1, hence the contour integrals § G m G m+p are computed exactly like those for F m , and the 
norm is unchanged. Note in particular that for q = 0, <p(x\q) oc x allows us to multiply 
F m (x) by any power of x. This gives us a lot of freedom in picking the model. 

Let us now present a simple example, with the same vector space E = IR <g> 1R <g> ... we 
have considered so far. 

4-2. Integrable models of discrete 2D Lorentzian gravity with polygonal tiles 

Let qi, qi, qk be generic distinct real numbers with \qi\ < 1. Let us introduce the 
function 



for m = 1,2,... We claim that this generates the components of a family of orthogonal 
vectors spanning E. The spanning property is clear by noticing that a triangular 

change of basis maps the canonical one K m (x) = x m onto the m 's& To see why they are 
orthogonal, we simply have to compute the following contour integral over the unit disc 

- ,Xx\q 1 ,...,q k )4> rn+p (-\q 1 ,...,q k ) = (-l) k f — x k \[ , (4.5) 



9 This is true because we assumed the q's to be generic, i.e. not related by any algebraic 
relation. It is easy to see for instance in the case k = 2, q± = q, qi = —q which is non-generic, 
that <f) m is an odd function of x, in which case only "half" of E is generated. 
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For p > 0, the integrand has no poles outside the unit disc, hence the integral vanishes, 
and the corresponding vectors are orthogonal. When p = 0, we find the normalizations 



/ 2m AX (l _ q jX )(x - qj) 
that reproduce the one-g result N(q) = a/1 — g 2 , and give for k = 2, 3 



(4.6) 



(l-g?)(l-g 2 2 )(l-g 1?2 ) 



1 + 9KZ2 



(1 - gg)(l - g 2 2 )(l - g|)(l - gig 2 )(l - gig 3 )(l ~ g 2 g 3 ) 

i + gig2 + gig3 + Q2Q3 - q^q^iqi + g2 + g3) - (<?i<?2g3) 2 



(4.7) 



The normalizations ( [4.6|) allow us to define orthonormal vectors i>( m ) = N(q\, qk)w^ m \ 
generated by F m (x\qi, g^) = N(qi, qk)(f>(%\qi, qk)- This produces an interesting 
candidate for the transfer matrix of Lorentzian gravity including triangles, squares,..., up 
to 2k + 2-gons, with the same fixed staircase boundary conditions as before: 

VIA) = — ^....^(^^(g,,...,^) (4g) 

Ui=i( 1 - ft^X 1 ~ QiV) ~ X Ui=i(Qi - x)(qi - y) 
To obtain a well-normalized generating function, we must take / of the form 

f (x) = aa-Kqi-qk) 2 ) , , 

^'-"^ j iV 2 (gi,...,g fe ) • (4 ' yj 

where yfg is the Boltzmann weight per boundary triangle, a free parameter of the theory, 
and the generating function for the transfer matrix elements finally reads 

= — k p- — (4.io) 

where 

w hJ (qu...,q k \X) = (-l)^ ^"/^-^" 3 " . ( 4 - n ) 

where the Oj = crj(qi, qk) are the symmetric functions of the g's defined by YliKjKki^j ~ 
x ) = ^2o<i<k(~ i y (J k-i x \ witn in particular a = 1 and a k = qxqz—qk- 
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The weights Wij Q4.ll ) are easily interpreted in terms of different Boltzmann weights 



for different types of polygons involved in the discrete world-sheet construction. Indeed, 
following the result of Sect 2.4 in the case of triangles and squares, the coefficients Wij are 
the Boltzmann weights of the following local configuration of i lower half-edges merging 
into j upper ones, all attached to the same point of the time-line: 



(4.12) 

with the obvious dual interpretation as an (i + j + 2)-gonal tile with two time-like edges 
and % lower- and j upper space-like ones. So the model whose transfer matrix is generated 
by ( |4.1Up is one of discrete 2D Lorentzian gravity involving polygonal tiles with 3,4,... up 



to 2k + 2 edges. 

Note that there would be a priori k{k + 3)/2 independent coefficients Wij in the most 
general symmetric generating function of the form ( [4.10Q , whereas here they are expressed 
in terms of only k + 1 independent parameters q±, qk, A. Hence, we are looking at a rather 
small-dimensional sub-manifold of the moduli space of 2D Lorentzian gravity, discretized 
with triangles, squares, ... up to (2/c + 2)-gons. The remarkable fact is that this manifold is 
integrable, as the corresponding transfer matrices commute with one-another. The physical 
meaning of the particular expression ( j4.11|) for the Wij weights is however not at all clear. 
In particular, for more that one qi, we cannot any longer map the model (as we did in 
Sect. 2.4 for one q) onto a "generalized curvature model" with only triangles and with 
additional curvature weights involving now nearest, next-nearest,... up to, say, /c'th-nearest 
neighbour interactions. Indeed, we can consider the most general model of triangles with a 
curvature weight a p for any succession of p up triangles (or p down triangles) and a weight 
g per triangle (and as before y/g for the boundary triangles). The integrable one-g case of 
Sect. 2 corresponded to the particular choice a p = a p_1 . Then is easy to show that the 
generating function 9(x,y\g, {a p }) generalizing ( |2.3| ) reads 

at I r u 1 A(gx)A(gy) 

^,v\9,M) = - gl _ A{llx)A(gy y (4.13) 

where A(u) = ^ p>1 a p u p is the generating function for the weights a p (for the one-q case, 
we simply had A(u) = u/(l — au)). The particular form ( [4.13| ) above then implies that 
the quantity 

gO{x,y\g, {a p }) 



1 +g6(x,y\g,{a p }) 
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A(gx)A(gy), (4.14) 



must be factorized into a function of x times the same function of y for some particular g. 
It is easy to check that for more than one qi, the generating function 9 qij .„ iQk (x,y\\) given 
by (|4 . 1 Q| ) does not satisfy this factorization requirement. 



Before we conclude this section, let us slightly extend the definition (|4.4|) to also include 
complex (and now non-generic) numbers q±,...,qk inside the unit disc, but such that <f> m 
remains real. A typical example would be to take q± = g, g2 = wq, <?3 = w 2 g,..., qk = u k ~ 1 q, 
where uj is a primitive /c-th root of unity. Then we simply have that <fi m (x\qi, qk) oc 
F m (x k \q k )/x k ~ 1 with F m as in (|2.6|). It is clear that the </> m 's only generate a fraction 
(roughly 1/k) of the infinite space E, but the content of the model is exactly the same as 
that of the one-g one. Another more interesting and possibly generic possibility would be 
to include both the g's and their complex conjugates in the list, namely take ( |4.4j ) with 
the 2k g's being gi, q%, g~i, g&. Note that the number of free real parameters in that 
case would still be 2k. We shall make use of this extension in the next section. 

Finally, we can extend the above construction so as to include other boundary condi- 
tions. The case of periodic boundary conditions is discussed in detail in Appendix A and 
corresponds to transfer matrices which are non- symmetric. In this case, the construction 
requires the use of two different (left and right) orthonormal bases. 

4-3. Continuum limits 

Let us now study the continuum limit of the new models (with k g's) defined in Sect. 
4.2. In this section, we are going to derive the rather negative result that no new scaling 
behaviour can be obtained from the models of Sect. 4.2, namely that the only well behaved 
continuum limit leads to the same result as in the one-g case of Sect. 2.6. 

To show this, let us first reexamine the case of one g. We start from the generating 
function flO) (WM iterated T times: 



«..t(», m = u.mr _ , (4.15) 



with 

1-Ag 2 
1-g 



UW=9-r^2- ( 4 - 16 ) 



To get a reasonable continuum limit, after taking T = r jet for some small parameter a, 
we need that both A = 1 + 0(a) and f q (X) = 1 + 0(a), while g tends to some critical value 
g c < 1 when a — > 0. From ( [4.16| ) above, this is only possible if g 2 = 1 + 0(a) too. So we 
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assume that q tends to 1 as q = e aV ^ = 1 — a\/X+0(a 2 ). Then from f q (X) = l + 0(a), 
we deduce at leading order in a that 

hm 9c{1 - X f ] = 1 A = l-2v / A^^« + 0(a 2 ). (4.17) 

o^o 1 - # c 

Once substituted back into the generating function 6 q (x,y\\) this gives 

lim d q (x, y\X) = . . 9cXy . r-^— , (4.18) 

«-o qy ' y| ; i-fi-^Xz + yHtl^c)^' V ; 

where we have identified the limiting Boltzmann weights ( |4.11| ) in ( [4. 1U| ) as 



w ,o -* 1, tw ,i = twi.o -> 1 - 9c wi,i ->• -(1 - 20 c ), (4.19) 

which remain all finite in this limit. 

Let us now use the same strategy to deal with the k-q model's continuum limit. The 
T times iterated generating function now reads 

e qi ,..., qk ; T & VIA) = (/ *-"" (A))T x y N \-^ , (4.20) 

rL=i(! - Qix)(i - Qiy) - a t Ui=i(Qi - x )(<ii - y) 

with 



as in ( |4.9|) and N(qi, g^) as in (|4.6|) . Setting again T = r/a and letting a — > (and 
g — > 0c) some critical activity per triangle), we must again have both f qi) ..., qk (A) = l + 0(a) 
and A = l+0(a). Let us moreover impose that when a — > 0, all the Boltzmann weights Wij 
of (|4.11|) remain finite. We must then take N 2 (qi, ...,qk) — > and therefore ((ft...?*;) 2 — > 1 



to be compatible with Eq. Q4.21| ) above. But since all the g's are real and less than 1 in 
absolute value, this implies the all qf — > 1 as a — > 0. Hence the natural generalization of 
<7 = e _aV ^ to the fc-g case reads 

q l = e ie - a % i = l,2,...,k, (4.22) 

for some small parameters ck, > and some signs = ±1. If we now do not impose any 
longer that the q^s are real, we can also choose pairs of conjugate complex numbers 

qi = uo l e~ ai , ft = uJie~ a \ (4.23) 

with c^j of modulus 1 and such that u 2 ^ 1. 
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Our choice is however strongly constrained if we insist in having finite Boltzman 
weights and it can be checked that, assuming that all the a^s are of the same order a, 
this implies that the norm N 2 (qi, q^) itself must be of order a. From the formula ( f4.7|) , 
we see for k = 2 and k = 3 that this in practice imposes that we can have at most one 
real qi tending to +1 (i.e. q = +1 in (|4.22|) above) and at most one real qi tending to —1 



(i.e. €i = —1). Finally, to reach a non trivial limit by setting as usual x = 1 — aX and 
y = 1 — aY, we need at least one qi tending to +1. We believe that these requirements are 
generic for arbitrary k. For k = 2, these constraints limit our choice to exactly two real 
gi's as in ( |4.22|) with e\ = +1 and e 2 = — 1. Writing 



ai = ay/A h i = l,2, (4.24) 

we may easily compute the scaling limit of ( f4.20|) in this case k = 2, with x = 1 — aX and 
y = 1 — aY, and we get 

C 

lim ad a , ao-T(x,y\X) = r-i— , , ; ; ; ^=r 5 (4-25) 

where = limA T , and C is a pre- factor independent of X and Y. This is simply propor- 
tional to the one-q result ( |2.29|) , and corresponds therefore to the same critical behaviour. 
As mentioned above, the choice 

qi = ue~ aVX , q 2 = ue~ aVX , (4.26) 

with ui a complex number of modulus 1 and such that uj 2 ^ 1, leads to a trivial limit. 
Indeed, in this case, iV 2 ~ |1 — uj 2 \ 2 VXa, and all the Boltzmann weights remain finite, 
but in the scaling limit where x = 1 — aX and y = 1 — aY, all dependence in X and Y 
disappears. In fact, as mentioned before, to retain a dependence on X and Y in the scaling 
limit, we must have at least one q tending to 1, which is not the case here. 

Applying now our constraints to the case k = 3, we need to choose exactly one of the 
qi's tending to +1 and the two others being complex conjugates, i.e. 

qi = e~ a \ q 2 =ue- a \ q 3 = ue~ a \ (4.27) 

with again \ui\ = 1 and uj 2 ^ 1. The normalization factor now becomes 

».2, \ 2\(l-uj)il-uj 2 )\ 2 a 1 a 2 , A N 

N 2 ( qi ,q 2 ,q 3 ) ~ - J \ , ; '_ ' \ 4.28 
2(a 1 + a 2 ) + (u + u)a 2 
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and we get finite limiting Boltzmann factors by taking 



«i = aVAi, i = l,2. (4.29) 
Now the scaling limit of ( |4.20| ) only retains the factors containing qi, and wipes the others 



out, so we finally get again the form (|4.25| ), with possibly different constants C and cf>. 

We are now ready to treat the general case. Depending on whether k is odd or even, 
we must take say 

k odd : qi = e" Q \ q 2i = Uie~ ai +\ q 2i+ i = u i e~ a *+ 1 
for % = 1,2, ...Ak- l)/2, 

(4.30) 

k even : q x = e Q \ q 2 = -e a2 , q 2i +i = w.e ai+2 , q 2i j r2 = Uj i e ai+2 
fori = 1,2, ...,(fc-2)/2, 

for some generic complex numbers cu^ of modulus 1, and such that ujf ^ 1. As only one 



q tends to 1, we are again led to a scaling limit of the form ([4.251) . Note that in the end 



of Sect. 4.2 we have already considered the case when qj = u : >~ 1 q, for j = 1, .., k and 
uj a primitive fc-th root of unity, and concluded then that the corresponding model was 
equivalent to that of k = 1, q± = q up to some minor rescalings (x — > x k , y — > which do 
not affect the scaling limit): this is in agreement with the present result, for a particular 
(non-generic) choice of the u^s as k-th roots of unity. 

To conclude, all our models have yielded the same and only non-trivial scaling limit 
for the two-loop correlator as the one-q one of Sect. 2.6. This is a manifestation of the 
rigidity of our integrability condition, in that all these integrable models share the same 
universality class. 



5. Conclusion 

We have revealed an interesting integrability structure underlying the simpler mod- 
els of two-dimensional Lorentzian quantum gravity and described how one can construct 
various extensions of these models while keeping them integrable. Among the models con- 
sidered so far, we have not found any for which the continuum properties of the Lorentzian 
geometries were different from those of the pure Lorentzian gravity case. We have concen- 
trated on a model of Lorentzian quantum gravity involving a higher curvature term which 
is equivalent to a model of Lorentzian quantum gravity interacting with a simple dimer 
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field. The fact that this model has the same continuum behaviour as pure Lorentzian 
quantum gravity is in contrast with the situation in Euclidean quantum gravity where the 



introduction of dimers changes the continuum behaviour of the geometrical system [11,19 



On the other hand our results are in compliance with the results of [[1511 , where no inter- 
action was seen when Lorentzian quantum gravity was coupled to an Ising spin system. 
However, one should bear in mind that when considered as a model describing matter 
coupled to Lorentzian gravity, our model contains only a subset of the possible matter 
configurations. It is our hope that our further investigations of the integrability structure 
revealed will enable us to solve exactly more realistic models of Lorentzian quantum grav- 
ity coupled to matter. In the light of the recent results of numerical simulations []T^| it 
would be particularly interesting to study the case of more than two Ising models coupled 
to Lorentzian quantum gravity. Interpreted as a model involving a higher curvature term 



our model is similar in spirit to the models considered in [|18| , where the effect of adding a 
higher curvature term to usual Euclidean quantum gravity was investigated. The result is 
the same for Euclidean as for Lorentzian triangulations. Adding a higher curvature term 
does not change the continuum physics for the geometrical system. 

We have furthermore proven the equivalence between Lorentzian triangulations and a 
certain type of random walks. This equivalence has allowed us to set up a dictionary con- 
necting concepts in Lorentzian quantum gravity to concepts in the theory of random walks. 
For instance, the loop-loop correlator of pure Lorentzian quantum gravity in the language 
of random walks became the generating function for large excursions. Furthermore, the 
integrability structure underlying the triangle-square model of Lorentzian quantum grav- 
ity could be understood using a simple block decomposition of the corresponding random 
walk. Finally, the equivalence between Lorentzian triangulations and random walks has 
provided us with an explanation why Lorentzian triangulations have fractal dimension two 
and why the continuum properties of the triangle-square model are the same as those of 
pure Lorentzian quantum gravity. It is very likely that a further investigation of the ran- 
dom walk picture will likewise tell us under which circumstances, if at all, we can expect 
new continuum behaviour to occur. We are also convinced that the above mentioned block 
decomposition will prove useful in revealing the integrability structure underlying the more 
realistic models of quantum gravity coupled to matter. 
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Appendix A. Non-symmetric transfer matrices and periodic boundary condi- 
tions 

It is instructive to derive the transfer matrix of the curvature model with periodic 
boundary conditions. Its features will allow us to enhance dramatically the construction 
of Sect. 4.2. The transfer matrix for the periodic model with curvature has a marked lower 
edge (that we will always represent as the leftmost one) , and always at least one upper edge, 
to avoid degeneration into the vacuum. The transfer matrix element Tfj r (g, a) between a 
row of i lower half-edges and j upper ones reads 



mi 777,2 



rrik 



E 



E 



k>l Sn T - i,n 1 ,..,n k >l 

Stn.f = j , m -L ,. . , m i > 1 , m ^, > 



ri2 



= E 



E 



gi+j oP{n r — 1)+E(m r — l)+2<5 mfe ,o 



fc>l Sn r —i,n 1 ,..,n, k >l 



k>l 



(A.1) 



where the boundary condition is taken into account by the fact that when the rightmost 
half-edge is a lower one, it receives a weight a as it is the neighbour of the leftmost one, 
instead of the a -1 given by the generic formula a mfe_1 at nik = 0. The generating function 
for ( |A.1| ) is readily found to be 



e^ r (x,y\g,a)= £ x^TfJ \g,a) 
i,3>l 



g 2 xy 



(A.2) 



(1 - gax)(l - ga(x + y) - g 2 (l - a 2 )xy)' 
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Let us perform the change of parameters g, a — > A, q of (|2.5|) and ( |2.19|) . The generating 
function ( |A.2| ) is transformed into 

g 2 (l - Xq 2 ) 2 xy 



O p q er (x,y\X) 



((1 - qx) - Xq(q - x))((l - qx)(l - qy) - X(q -x)(q- y)) 
A ^ (1_9) (l-qx)r+° X (1 - ^ 



r,s>l 



(A.3) 



A £ P m (x| ? )A m - 1 Q m (y| ? ). 



m>l 



where 



P m (x|g) = v 7 !-? 2 



x(q — x) 



m— 1 



(1 - qx) m+1 



*.«.)- d-rti E ^'^-^ <<--(££)">■ (A ' 4 ' 



7", S > 1 

r+ s = m. + 1 



We know already that the periodic case with curvature is solvable, as we have computed 
its two-loop correlator above. The reason why is extremely simple: the functions 

P m (x\q) = J2* l 4 m \ 

^ , , (A-5) 

Q m (y\q) = J2y jw i > 
i>i 

generate the components of two different left and right vectors w^ m \ n/ 771 ) for T per , that 
turn out to be orthonormal to one another, i.e. w;( m ) • w^ m ) = S m m >. It is indeed easy 
to show that the contour integral over the unit circle: § Pm^x^Qmi (1 / x\q)dx / {2inx) = 
dm,m'i by use of Cauchy's residue formula. 

Note however that these are not eigenvectors of T per but we have the "diagonal" 
mapping 

T per w {m) = X m w {m \ (A.6) 

This translates immediately into the following transfer matrix relation (as before we define 
TP er (X) = TP £r (g, a) with g,a—>X, q) 

TP er (X)TP er (X') =TP er (XX'). (A.7) 

Note that this relation is also satisfied by T q (X) ® T q (X), thanks to (|2.23| ), and justifies a 
posteriori our previous calculation of the two-loop periodic correlator (p = 2 in ( |2.40| )). 
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Note finally that the two sets of (left and right) vectors span the same space E, as Q m {y) oc 
y for small y, and there is a triangular change of basis from the -F m 's to the P m 's or Q m 's. 

Let us draw some general conclusion from this example. A score of other models can 
be constructed by means of two different sets of left and right vectors, provided those are 
mutually orthonormal. Given any such pair generated by say P m (x) and Q m (y), we get a 
set of transfer matrices T(A) with spectral parameter A, generated by 

9(x,y\X) = T h3 {X)x^ = /(A) Pm{^ m - X Qm{yl (A.8) 

i,j>l m>l 

for some arbitrarily chosen function /(A). The vector spaces spanned by the left and right 
vectors need not be the same, we simply need the vectors to be mutually orthogonal: 
§ dx/(2mx)P m (x)Q m i(l/x) = 8 m>m i. This implies a multiplicativity relation 

T(A)T(A') = /( / A ( ) A / A ( / ^ ) T(AA / ). (A.9) 

which makes the model trivially integrable. 

In the light of this, it is easy to find a non-symmetric generalization of the k g's model 
of Sect. 4.2. Taking 

k ' \ra- 1 



(Qj ~ x ) 



p m (x) M k ./• ] T — 

(A. 10) 

Qm(y) = {-i) k - x M k (( qiq2 ... qk r - II (TT^r) ■ 

as generating functions of the left and right vectors respectively, it is easy to prove that 
these are mutually orthonormal, provided we take 

1 , .vt f dx TT 1 



Ml J 2i7rfJ i (l-q j x)(q j -x) 

k k , k , (A- 11) 

~[ t\ 1 ~ qiq * i=T * " Qi 

by Cauchy's formula. For k = 1, this gives Mi = 1 — g 2 , in agreement with ([A.4j ), while 
for k — 2, 3 we have 

M 2_ (l-gl)(l-g 2 2 )(l-glg2) 
<?l+<?2 

M| = (A. 12) 

(1 - g 2 )(l - gj){l - g§)(l - gig 2 )(l - g 2 g 3 )(l - gigs) 

g 2 + g| + gf + gig 2 + g 2 g3 + gig3 - gig 2 g3(gi + g 2 + g3) - q\q\ - q\q\ - q\q\ ' 
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The vectors ( [A .10 ) lead to the transfer matrix generated by ( |A.8| ) 



9(x, 7/| A) = xg(X) J — 



n, =1 (l - q,x) -A[[ j=1 qMj - x) (A. 13) 

nj=i(! - - - a nj=i(^' - %)(qj - y) 

where g(X) is an arbitrary function of A (in which we have absorbed various normalization 
factors). We suspect the transfer matrix of ( |A.13| ) to be the periodic boundary condition 



version of (|4.8| ) ( fO| ). Note that here the vectors generated by ( |A.10|) both span the same 
space E (they are proportional to x and y respectively). 

Let us finally mention the following generalization of ( |A.10|) , based on the general 
property that if P m (x) and Qm(y) are generating functions of mutually orthonormal 
vectors, then for any function R(x), the vectors generated by </> m (x) = R(x)P m (x) and 
i^miy) = Qm(y)/R(^/y) are also mutually orthonormal (we must of course assume that 
both R(x) and 1/ R(l/y) have convergent series expansions). We may apply this recipe to 
( |A.10| ) with the function R(x) = 1 / (0^=1 (1 ~ 9.7 ^O^ -1 ^ f° r some given integers pj > 1. 
But l/R(l/y) = (—l) p ~ k n^Life' —y) Pj ~ 1 /y p ~ k i where we have defined p = J^Pji has no 
good series expansion around y = 0. In order to fix this and arrange for both the left and 
right vectors to generate the same space E as above, we are led to take 



N m — 1 



<j> m (x) = s k x TT } qj x s 

rmv ) K 11 j_ X )m+ Pj 

i=i J 



(A.14) 



The coefficients a\ are fixed by requiring that i^m(y) = 0(y), and Sk is the normalization 
factor ensuring orthonormality, with the result Sk = defined in ( |A.11| ). Note that for 
any polynomial P, P(l/y) is automatically orthogonal to <p m (x), as no pole lies inside the 
unit disc in the corresponding Cauchy integral § 4> rn (x)P(x)dx/ (2iiTx) = for all to > 1 
(the x in the denominator is cancelled by the one in factor of </> m ). 
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